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M* H. Vavra 



ABSTRACT 



The paper describes a graphical method for the 
solution of the matching problem in gas turbine plants. 
The method is used to investigate the conditions in a 
closed-cycle system similar to that of the ML-1 nuclear 
gas turbine plant of Aerojet-General Nucleonics. It is 
shown that proper matching of the turbomachines with a 
given system is of great importance to obtain the best 
possible performance. 
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NOMENCLATURE 
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KW 
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ATM 

V 
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e 

M* 

5 

p 



constant for pressure drop calculations 
mass flow rate (Ibm/sec) 

total pressure (psia) 

760 mm Hg = 14.7 psia 

pressure ratio 
speed (rpm) 
power (KW) 

total absolute temperature 
288°K = 518.4° R 



velocity 

P ^ P ATM 



T/T 



ATM 



leakage flow coefficient 
pressure drop coefficient 
mass density of fluid 



Subscripts 



1. refers to compressor inlet 

2 M 11 compressor discharge 

3 11 " turbine inlet 

4 M 11 turbine discharge 

C " 11 compressor 

T " 11 turbine 

5 n n shaft power 

HP 11 n high pressure side of loop 

LP n " low pressure side of loop 
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INTRODUCTION 



In closed-cycle gas turbine plants it is quite difficult to evaluate 
the influence of the different components on the performance of the whole 
system. These difficulties arise because of the shifting of the operating 
points of the turbomachines, that produces changes in pressure ratios and 
mass flow rates. Without a computer it is necessary to perform lengthy 
calculations with methods of successive approximations to determine these 
changes. Even if a computer program is available the graphical method for 
the solving of the matching problem, which is discussed in the following, 
has merits since it gives a good insight how the performance of the dif- 
ferent components affect the operating conditions of the plant. 

The method will be applied to the simple cycle of the nuclear ML-1 
plant of Aerojet-General Nucleonics of Fig. 1. It is assumed that the 
performance maps of the turbomachines are given, that the pressure drops 
in the system be known, and that the plant operates between given temp- 
erature levels. More complicated cycles with multiple compressors and 
turbines, having intercooling and reheat, can be investigated with the 
same general principles that are outlined in this paper. 

If m denotes the mass flow rate entering the compressor of Fig. 1, 
the turbine flow rate m^ is 

m T = (1 - |i) m c (lbm/sec) (1) 

where |j, represents that fraction of the flow rate which escapes from 
the high pressure to the low pressure side of the loop, either through 
the bypass control valve, or through seals and leakage paths. The fluid 
that escapes from the loop to the surroundings is replenished by a make- 
up system, controlled by the pressure of the cycle. The power required 
for the compression of the make-up gas will be ignored. 

The pressure p^ ahead of the turbine can be expressed by 

P 3 - (1 - § Hp ) P 2 (2) 

where p^ is the compressor discharge pressure and a factor that 

depends on the pressure drops in the whole high pressure loop. 

The turbine discharge pressure p^ can be expressed similarly by 

P 4 “ (1 + W P 1 (3) 
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with the compressor inlet pressure p^, and the pressure loss coefficient 
of the low pressure loop of the plant. The pressure ratio r^ * p^/p^ 
of the turbine is then related to the pressure ratio r^ * $ 2^1 t * ie 
compressor by 

_ ^ ~ ... 
r T ' (1 + ^ Lp ) r C ( } 

Since the turbomachines of the plant under discussion are directly 

coupled, compressor speed N p must equal the turbine speed N . This con- 

L* 1 

dition, and Eqs. 1 to 4, must be satisfied, and made compatible with the 
performance of the turbomachines, as given by their maps, to solve the 
matching problem of the system. 
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PERFORMANCE MAPS OF TURBOMACHINES 



The performance of a compressor that operates with a particular fluid 
can be represented by 



and 







KW C 
6 - 1 / 0 ! 



( a c / 6i 

2 V 5„ 



/e n 






(5) 



( 6 ) 



where KW represents the power in KW to drive the machine. The quantities 
0 and 5 are defined by 



0 



T/T 



ATM 



5 



P/P 



ATM 



The arbitrarily chosen reference values T Am> , ■ 288°K a 518.4° R, and 

ATM 

P atm = 760 mm Hg B 14.7 psia, make it possible to express the referred flow 
rates, speeds, and power requirements with conventional units. The subscripts 
of 4 and 8 correspond to the stations in Fig. 1 where the total absolute 
temperatures T and the total absolute pressures p occur. 

In contrast to usual practice the referred flow rate of the compressor 
is taken with respect to the discharge pressure and the inlet temperature, 
for convenience in the later calculations. If the flow rate 
given instead of the quantity used in Eqs. 5 and 6 , the latter can be es- 
tablished also, since at a particular referred speed N^//©^ and a referred 
flow rate 3 ^ / 5 there occurs only one pressure ratio P 2 /P^ # Then 

"c /e i . *C / 9l fl 
S 2 6 ^ P 2 



It is of interest to note that the referred power of Eq. 6 is the 
power necessary to drive a compressor whose inlet conditions correspond 
to the standard atmosphere defined above, and whose speed is such that 
the same velocity triangles occur as in the actual machine. 
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For the present study it is advantageous to express the performance 
of a turbine by 



/9o Po N 

— - f 3 < 7 > 

and 

kw t _ P 3 N t 

sy^ - ■ ^ T h ’ 7e7 ) (8) 

The referred turbine flow rate of Eq. 7 and the referred shaft power 
of Eq. 9 represent the flow rate and power, respectively, that the turbine 
would have if the inlet conditions of the gas were those of the standard 
atmosphere. 

Two turbocompressor sets were developed for the ML-1 plant. One, 
which henceforth will be called TC-set M a", has a two-stage centrifugal 
compressor and a two-stage axial turbine, operating at a nominal speed 
of about 18,000 rpm at 118°F 3 48°C compressor inlet temperature, and 
1200°F = 650°0* turbine inlet temperature. The other set ("B") has an 
eleven-stage axial compressor and a two-stage axial turbine with a de- 
sign speed of about 22,000 rpm at the above-listed temperatures. Both 
sets operate with nitrogen as working fluid, to which 57 0 oxygen is added 
to prevent nitriding of the reactor materials. 

Although the method of this paper can be used to obtain matching 
solutions at arbitary speeds and cycle temperatures, the present investi- 
gations will be restricted to the above-mentioned temperature levels and 
the design speeds of the TC sets. Then, for TC set "A 11 : 



\ 
/ 9 



N, 



- 77 - - 17,170 rpm, and, 



T _ 



/e, 



= 10,130 rpm 



For TC set "B": 



N c 

/e 



* 7 — - 20,620 rpm, and, 



/e, 



12,160 rpm 



For these referred speeds the functions f^ of Eq. 5 of the compressors 
of both TC sets are presented in Fig. 2. The functions f^ of Eq. 7 for 
the turbines are shown in Fig. 3. It can be noted that both figures have 
the same double logarithmic scales since this representation is required 
for the present method. Fig. 4 and 5 give the powers required by the 

compressors, and produced by the turbines, respectively, at the referred 
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design speeds, in accordance with Eqs. 6 and 8. The data in Figs. 2 
and 5 pertain to the state of the TC sets as they were initially de- 
signed and fabricated, and represent the best estimate of their initial 
performance obtained by a number of tests with different instrumenta- 
tion at various conditions. 



SYSTEM PRESSURE LOSSES 

Pressure drops due to friction, flow separations, or mixing can 
be expressed by 



Ap = f p V 
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where the dimensionless coefficient f may vary with Reynolds number. 
The quantities p and V are the average mass density of the fluid 
and the average flow velocity in the element where the total pressure 
is reduced by Ap. Since V =* m/( p A ), where A is the cross-sectional 
flow area, there is with p oc p/T, 



A£ 

P 




( 10 ) 



The quantity c is a constant for a particular element of a system. 
The total pressure drop in the low pressure side of the loop of Fig. 1 
may be expressed by using Eqs. 3 and 10 



p 4 - Pi 



^LP Z C 



.2 

m T. 

l 

i 2 

P, 



where the index i refers to the different elements that produce pressure 



drops. For small values of £ there is p^ 



'l’ 



and with m = m. 



- , m~/9 



T. 



s LP -^i; ies j 



an 



ATM 



If the temperatures in the loop remain nearly constant, Eq. 11 may be 
written as 



^LP 




m c /e 




( 12 ) 



where k will be a constant also. 
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With a similar approach the pressure drop coefficient £ of 
Eq. 2 for the high pressure loop can be expressed by 

2 



=HP 



k HP U - \J 



( 13 ) 



if it is assumed that the leakage factor p, of Eq. 1 is entirely due to 
flow through the bypass control valve. 

The measured pressure drop coefficients F and F of the ML-1 

^LP nr 

loop are plotted in Fig. 6. With these data the coefficients k^ p and 
k of Eqs. 12 and 13 were determined, and are shown in Fig. 6 also. 

nr 

It can be noted that k is very nearly constant, whereas k decreases 
somewhat with increasing referred flow rates. This condition is prob- 
ably due to the widely changing temperatures in the high pressure loop. 
The results of Fig. 6 will be used for the subsequent system analysis. 



GRAPHICAL SOLUTION OF MATCHING PROBLEM 



Eq. 1 will be expressed with the referred flow rates of the turbine 
and the compressor, as used in Eqs. 5 and 7> 



With Eq« 2 



m T ^3 _ m c p 2 ,, v x 3 
~ 3 



V e 3 V 9 1 1 



S 2 1 



^ h 

^HP / T 1 
V 



(14) 



Further, from Eq. 4 



r T - r c 1 ' ? HP 



1 + 



(15) 



■LP 



By taking the logarithms of the terms of Eqs. 14 and 15, and by 
considering the values thus obtained from Eq. 14 as vector components in 
the directions of the unit vectors i, those from Eq. 15 as vector com- 
ponents in the directions of the unit vectors j of Figs. 2 and 3, and 
performing a vector addition, there is 
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+ j i°g r c 



( 16 ) 



With this procedure it has been possible to separate completely 
the turbine performance data (two terms on the left-hand side of Eq. 16) 
and the compressor performance data (first two terms on the right-hand 
side of Eq. 16) from those parameters which are affected only by the 
conditions in the system that connects the turbomachines, or by the temp- 
eratures at which heat is added to or dissipated from the cycle. The last 
two terms of Eq. 16 will therefore be called the system transfer vector S, 
whereas the two terms on the left-hand side are denoted by T, and the first 
two terms on the right-hand side by C. Then, the vector relation 



satisfies the matching condition of the turbine and the compressor for a 
particular loop. Since an arbitrary constant vector can be added to T 
and C without invalidating Eq. 17, the system transfer vector S indicates 
simply how points with the same coordinates in the turbine and compressor 
maps of Figs. 2 and 3 must be displaced in i and j directions to account 
for the systems characteristics in the matching solution. Possible solu- 
tions are then obtained for = N^, > or 



where the ratio T^/T^ is again depending only on the imposed temperature 
limits, and not on the performance maps. Hence, if performance maps of 
the turbine and the compressor that are drawn on transparent paper are 
displaced in accordance with the vector S, the possible operating points 
of the machines can be found at the intersection of the characteristics 
that satisfy Eq. 18. 

From Eq. 16 



T = C + S 



(17) 




(18) 




+ log (1 - n) - j 



log (l+| Lp ) 



(19) 



+ (j - i) log (1 - | Hp ) 
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Fig. 7 shows a diagram for the determination of the vector S. With 

the unit vectors i and j in the same directions as the corresponding vectors 

in Figs. 2 and 3, the vectors S obtained from Fig. 7 have their origins at 

a particular station, say X^ , in the turbine flow diagram and their termini 

at the station X in the compressor flow diagram which has the same coordinates 
L* 

as X T (see Fig. 8) . 

From Fig. 7 it is evident how the different system parameters influence 
the vector S. For instance, leakage or bypass flow has the same effect as 
a reduction in turbine inlet temperature. If T^/T^ is the actual ratio of 
turbine and compressor inlet temperatures, a certain bypass flow, expressed 
by the factor M-, is effectively reducing this ratio to 




(1 - ^) 2 



Hence, for * 1200°F and * 118°F a leakage flow rate of 5 % of m 
has the same effect as a reduction of the turbine inlet temperature from 
1200°F to about 1040°F. This example shows that output control by means 
of bypass flow is very effective, but it shows also the detrimental effect 
of internal leakage flows from the high pressure to the low pressure side 
of the loop. 



EXAMPLE 

In Fig. 8 are shown the compressor and turbine flow maps of Figs. 2 
and 3 displaced in accordance with the system transfer vector S from a 
station X^ to a station X^, , which have the same coordinates in the respec- 
tive maps. This shift of the coordinates establishes the graphical solu- 
tion of Eq. 17, Since the pressure drop coefficients F and F depend 

HP °LP 

on the flow rate and pressure ratio to be determined, the actual vector S 
must be obtained with a method of successive approximations. The pressure 
drop coefficients were determined with Eqs. 12 and 13, using the average 
values k^p = 0.0036 and k^ p = 0.0525 from Fig. 6. The matching is carried 
out for T^ = 1200°F, T^ - 118°F 5 and a bypass flow of 37, of the compressor 
flow rate, at the nominal design speeds of the TC sets. Table 1 shows the 
match point conditions obtained from the points of intersection of correspond- 
ing characteristics in Fig. 8. 
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Table 1 



Match Point Solutions from Fig. 8 



TC Set 


A 


B 


rn /e, 

~ (lbm/sec) 

2 


1.362 


i 

? 

i 

i 

1.375 > 

i 

( 


r C ' p 2 /p l 


2.725 


2.285 I 

] 


m_/e, 

g (lbm/sec) 

°3 


2.460 


2.480 


r T = P 3 /P 4 


2.38 


2.00 


? LP 


0.049 


0.036 


^HP 

! 


0.090 

1 

L 


0.091 



With these data it is possible to determine the shaft powers of the 

TC sets from Figs. 4 and 5, where the match points of Table 1 are denoted 

by M a and M_. . 

J A B 

The power of the TC sets will be calculated for a compressor inlet 
pressure p^ = 117.6 psia or 5^ 5:0 8, since higher pressures p^ produce 
excessive stresses in the reactor vessel. The results are listed in 
Table 2. 
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Table 2. Power of TC Sets for Match Point Solutions of Fig. 8. 



RW C 



(KW), (Fig. 4) 



218.0 



149.6 



KW„ 



& 3 /e 3 



(KW), (Fig. 5) 



60.2 



48.0 



& 3 = (i - F Hp ) 6 2 - (1 - F Rp ) &1 - 21.80 18.28 



KW„ 



KW. 



(KW) 

(KW) 



2135 



1842 



1429 



1264 



KW S - KW t - KW C 



(KW) 



293 



165 



DISCUSSION 

The specifications for the TC sets called for a shaft power of 
about 480 KW at the design point. The reasons for the power deficiency 
of the sets become evident from Figs. 2 and 3 where the match point solu 

tions of Table 1 are marked by M, and M for the TC sets A and B. 

A B 

At the specified design point the referred turbine flow rate should 
be 2,33 lbm/sec. Fig. 3 and Table 1 show that the actual flow rates of 
the TC sets A and B exceed this value by about 5%% and 6%7 0 , respectively 
at the match points. The compressor of TC set A is therefore forced to 
operate at an efficiency of about 807, whereas it could have about 847, if 
the flow rate were lower. Although the efficiency of the compressor of 
TC set B is generally lower than the design estimates, a reduced turbine 
flow rate would make it possible to operate the compressor at an 
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efficiency of about 78.5% instead of at about 76. 5% at the match point. 
The increased flow rates also cause higher pressure drops on the high 
pressure side of the loop since £ is now about 9% instead of about 67o 
in accordance with the design point data. 

The turbines of the TC sets have flow areas of about 10 sq.in. at 
the discharge of the rows of blades, where the throat widths of the flow 
channels between neighboring blades are between 0.15 and 0.20 in.. If 
the blade profiles have machining tolerances of + 0.002 in., changes 
in flow area of between 2 and 47, are possible within the specified tol- 
erances. The writer is also of the opinion that the usual calculating 
methods for turbines must not be applied to small turbines without giv- 
ing special attention to the dimensioning of the flow areas. It is 
common practice to use the same loss coefficients for the efficiency and 
the flow area calculations. However, since the flow rates will be de- 
pending on the flow conditions at the throat of the blades, rather than 
on the conditions after them, it is necessary to apply smaller losses 
for the area calculations. If this is not done the flow areas will be 
larger than necessary. To investigate these conditions a computer pro- 
gram was established Til*, which uses the blading losses of T2]. Only 
the so-called profile losses are applied to the flow area calculations 
by establishing an equivalent displacement thickness of the boundary 
layers at the throat areas of the blades. This program, which also takes 
account of after-expansions if the flow is chocked at the discharge of a 
blade row by applying the relations of T 3 ~1 , iterates accurately the flow 
conditions after each row of blades, and produces complete performance 
maps for different speeds and arbitrary fluids. Although at present re- 
stricted to two-stage turbines, designs with an arbitrary number of 
stages can be investigated with slight modifications of the program. 

Another computer program P4] was set up to investigate means for 
improving the performance of the axial compressor of TC set B, and of 



Numbers in brackets designate References at end of paper. 
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shifting its characteristics to higher flow rates by increasing the 
stagger angles of the guide vanes. The program is based on the cascade 
data of T5] and takes account of the tip clearance, end, and secondary 
losses with the relations of f6]. The losses in the entrance and dis- 
charge passages of the compressor are considered also, and the program 
calculates the performance of each row of blades, and iterates exactly 
the flow conditions after each row. The deterioration of the velocity 
profiles toward the compressor discharge is taken into account by in- 
troducing blockage factors due to wall boundary layers, and work-done 
factors because of pointed velocity profiles. The two last-mentioned fa 
factors were obtained for the machine at hand by applying a variational 
process until the calculated performance matched the available test re- 
sults. Additional work in this direction is presently carried out at the 
Propulsion Laboratory of the United States Naval Postgraduate School, 
where tests are being made with a completely instrumented twelve-stage 
axial compressor. 

In the following it will be demonstrated what improvements of the 
TC sets A and B are possible with the existing compressors, by simply 
reducing the turbine flow rates to levels where the compressors can 
operate near their maximum efficiency. The chosen operating points are 
denoted by 0^ and 0 g in Fig. 2. The modified turbines are supposed to 
have about the same efficiency as the original designs, namely, 817o for 
TC set A and 807o for TC set B. The results are shown in Table 3 for S^=* 8 

Hence by simply reducing the turbine flow areas, the shaft powers 
can be increased by about 257o and 657, for TC sets A and B, respectively. 

If at the same time the turbine efficiency of both sets could be im- 
proved to 857o the shaft powers would be 455 KW and 383 KW, respectively, 
for sets A and B. 

These results clearly show that the proper matching of turbines 
and compressors in closed-cycle gas turbine plants is of utmost importance 
Even machines with high efficiencies will not produce the required plant 
power if a mismatch prevents them to operate at the optimum efficiency. 
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Table 3 Operating Performance of TC Sets for Turbines with Reduced 
Flow Areas 



TC Set 


A 


B 


/e, 

~ (lbm/sec) 

°2 


1 

1.15 


1.15 


r C = P 2 /p l 


J 

2.900 


2.665 


^HP 


0.065 


0.065 


^LP 


0.040 


0.038 


r T = P 3 /P 4 ( Ec l- 


2.605 


2.410 


KW (KW) 


1671 


1439 


KW X (KW) 


2036 


1714 


KW g » KW t - KW C 


365 


275 
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M. Ho VAVRA Fig* 5 Referred Power of Turbines of Turbocompressor Sets. 
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